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Abstract. We consider random g-state Potts models for 3 < g < 8 on the square lattice where the fer- 
romagnetic couplings take two values Ji > J2 with equal probabilities. For any q the model exhibits a 
continuous phase transition both in the bulk and at the boundary. Using Monte Carlo techniques the 
surface and the bulk magnetizations are studied close to the critical temperature and the critical expo- 
nents /3i and /3 are determined. In the strip-like geometry the critical magnetization profile is investigated 
with free-fixed spin boundary condition and the characteristic scaling dimension, Pi/v, is calculated from 
conformal field theory. The critical exponents and scaling dimensions are found monotonously increasing 
with q. Anomalous dimensions of the relevant scaling fields are estimated and the multifractal behaviour 
at criticality is also analyzed. 
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1 Introduction 

The presence of quenched i.e. time independent disor- 
der could modify the cooperative behaviour of physical 
systems with many degrees of freedom. In classical sys- 
tems, where thermal fluctuations dominate quantum fluc- 
tuations the effect of disorder in the pure system phase- 
transition point can be analyzed by relevance-irrelevance 
criterions. For second-order transitions, according to the 
well known Harris criterion [1], disorder appears to be a 
relevant perturbation which moves the random system to- 
wards a new fixed point when the specific heat exponent a 
\ of the pure system is positive. In the other situation, a < 
0, the disordered system remains in the pure model uni- 
versality class. The two-dimensional random-bond Ising 
model (RBIM) corresponds to the marginal case. It has 
been extensively studied in the 80's (for reviews of the- 
oretical and numerical studies, see Refs. [2] and [3], re- 
spectively) . The effect of randomness on first-order phase 
transitions was considered later. Imry and Wortis argued 
that quenched disorder could induce a second-order phase 
transition [4]. This argument was then rigorously proved 
by Aizenman and Wehr, and Hui and Berker [5,6]: In two 
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dimensions, even an infinitesimal amount of quenched im- 
purities changes the transition into a continuous one. 

The random bond Potts model (RBPM) is the para- 
digm of systems the pure version of which undergoes a 
second-order or a first-order transition, depending on the 
number of states, q, per spin [7]. In two dimensions, the 
second-order regime (7 < 4 has been considered by a num- 
ber of authors, using perturbative field-theoretical tech- 
niques [8,9,10,11,12,13] or Monte Carlo (MC) simulations [14,15,16]. 
On the other hand, in the first-order regime, g > 4, confor- 
mal perturbation techniques can not be used around the 
pure model transition point and the resort to numerical 
calculations becomes essential. Both Monte Carlo simula- 
tions and Transfer Matrix (TM) techniques, combined to 
standard Finite Size Scaling (FSS) [17,18,19,20,21,22] and 
conformal methods [23,24,25,26,27] were used at the ran- 
dom fixed point of self-dual disordered models to study 
bulk critical properties. 

The surface properties of dilute or random-bond mag- 
netic systems were on the other hand paid less atten- 
tion. Generally surface quantities, such as magnetization, 
energy-density, etc. are characterized by a different set 
of scaling dimensions, than their bulk counterparts. For 
example in the pure Ising model, bulk magnetization van- 
ishes as m ~ t^, with f3 = 1/8, whereas for the surface 
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magnetization the decay-law, toi ^ t^'^ , involves the sur- 
face exponent f3i ~ 1/2, where t denotes the reduced 
temperature. Quite generally, the scaling laws involving 
surface and/or bulk exponents can be deduced from the 
assumption of scale invariance. For example the singular 
part of the bulk, /b, and surface, /s, free-energy densities 
in a d-dimensional system behave under a scaling trans- 
formation, when lengths are rescaled by a factor 6 > 1, 
/' — l/b, as 

Mt,h) = b-''MbyH,by'^h), (i) 

Mt,h,hs) = b-^^-^^MbyH,by^h,by'^^hs). (2) 

The whole set of bulk and surface critical exponents can 
be expressed in terms of the anomalous dimensions yi as- 
sociated to the relevant scaling fields [28] (temperature 
t, bulk h and boundary hs magnetic fields), for example 
P= {d- yh)/yt and (3i = {d - 1 - yhMvt- 

The (1,1) surface of the disordered Ising model on a 
square lattice has recently been investigated through MC 
simulations by Selke et al. [29,30]. (For a related study 
of the critical behaviour at an internal defect line in the 
disordered Ising model, see Ref. [31].) The critical expo- 
nent Pi was found robust against dilution keeping the pure 
system value /3i = 1/2 and no logarithmic correction has 
been observed, in contradistinction with the correspond- 
ing bulk behaviour. The surface properties of the 8-state 
RBPM were also considered in Refs. [19,27]. 

In this paper, we report extensive MC and Transfer 
Matrix studies of the critical behaviour of both the sur- 
face and bulk magnetizations of the disordered Potts fer- 
romagnets for different values of 3 < g < 8. Our study ex- 
tends previous investigations in several directions. First, 
we investigate the temperature dependence of the bulk 
and surface magnetizations and calculate the critical ex- 
ponents (3 and /3i. Second, we consider strip- like systems 
with fixed spin boundary conditions (BC), determine the 
magnetization profile at the critical temperature and cal- 
culate the scaling dimensions Xb = 2 — y^ = P/v and 
Xi = 1 — yh^ = Pi/v from predictions of conformal in- 
variance. Our third investigation concerns the possible 
multifractal behaviour of the correlation function and the 
critical magnetization profile. The n-th moments of both 
quantities are found to follow predictions of conformal in- 
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are ob- 



variance and the scaling dimensions x\ 
tained n-dependent. 

The structure of the paper is the following. In Sec- 
tion 2, we present briefly the model and the simulation 
techniques. Section 3 is devoted to the approach to crit- 
icality, while in Section 4, magnetization profiles in the 
transverse direction of strips with fixed- free boundary con- 
ditions are computed. Multifractality is studied in Sec- 
tion 5 and a discussion of the results is given in Section 6. 



2 Model and algorithms 

2.1 The random-bond Potts model 

We consider Potts-spin variables, ai^k G 1, 2, . . . , g on the 
sites of a square lattice with I — 1, 2, . . . , L columns and 



k — 1,2, . . . , K rows, with independent random nearest- 
neighbour ferromagnetic interactions Jik and J/j, in the 
horizontal and vertical directions, respectively, which have 
the same distribution and could take two values, Ji > J2, 
with equal probabilities: 

T^iJik) = \s{Jik - Ji) + ^S{Jik - J2). (3) 



The Hamiltonian of the model is thus written 



n^ 



l,k 



{JikS. 



<^l,k:(^l + l,k 



^lk^<7ik,'yi,k-\-l) 



(4) 



In the thermodynamic limit L, K 
self-dual and the self-duality point 



00 the model is 



[exp(Ji/fcsT,) - l][exp(J2/fcsrc) - 1] - g 



(5) 



corresponds to the critical point of the model if only one 
phase transition takes place in the system. This assump- 
tion is strongly supported by numerical calculations. 

The degree of dilution in the system can be varied 
by changing the ratio of the strong and weak couplings, 
r — J1/J2. At r = 1, one recovers the perfect g-state 
Potts model, whereas for r ^ 00 we are in the percolation 
limit, where Tc = 0. The intermediate regime of dilution 
1 < r < 00 is expected to be controlled by the random 
fixed-point located at some r = r*(q). 



2.2 Monte Carlo simulations 

For the simulation of spin systems, standard Metropo- 
lis algorithms based on local updates of single spins suf- 
fer from the well known critical slowing down. As the 
second-order phase transition is approached, the correla- 
tion length becomes longer and the system contains larger 
and larger clusters in which all the spins are in the same 
state. Statistically independent configurations can be ob- 
tained by local iteration rules only after a long dynamical 
evolution which needs a huge number of MC steps and 
makes this type of algorithm inefficient close to a critical 
point. 

Since disorder changes the transition of the Potts model 
into a second-order one, the resort to cluster update al- 
gorithms is more convenient [32,33]. These algorithms are 
based on the Fortuin-Kasteleyn representation [34] where 
bond variables are introduced. In the Swendsen-Wang al- 
gorithm [35] , a cluster update sweep consists of three steps: 
Depending on the nearest neighbour exchange interac- 
tions, assign values to the bond variables, then identify 
clusters of spins connected by active bonds, and eventu- 
ally assign a random value to all the spins in a given clus- 
ter. The Wolff algorithm [36] is a simpler variant in which 
only a single cluster is flipped at a time. A spin is ran- 
domly chosen, then the cluster connected with this spin is 
constructed and all the spins in the cluster are updated. 

Both algorithms considerably improve the efficiency 
close to the critical point and their performances are com- 
parable in two dimensions, so in principle one can equally 
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choose either one of them. Nevertheless, when one uses 
particular boundary conditions, with fixed spins along some 
surface for example, the Wolff algorithm is less efficient, 
since close to criticality the unique cluster will often reach 
the boundary and no update is made in this case. In the 
following, we will consider two different series of simula- 
tions, one with free BC where the Wolff algorithm will be 
used, and the other with fixed-free BC for which we have 
chosen the Swendsen-Wang algorithm. 



3 Approach to criticality 

In this Section we consider square shaped systems, where 
L and K are equal, with L ranging from 40 to 640 to 
check finite size effects. In the vertical direction we im- 
pose periodic boundary conditions, thus (Ji^k+i = <7i,i, 
for I = 1, 2, . . . , L, whereas in the horizontal direction the 
boundary spins ai I = 1 and I — L are free. Thus we have 
a pair of (0, 1) surfaces, obtained by cutting bonds along 
the vertical axis of the system. 

According to numerical studies about the bulk quan- 
tities of the random model the finite size corrections are 
very strong unless the calculations are performed close to 
the random fixed-point [27]. The approximate position of 
r*{q) is listed in Table 3 for different values of q, as ob- 
tained from the maximum condition of the central charge 
of the model [27,37]. Our simulations were performed at 
these fixed-point values of the dilution, but for comparison 
we have also considered systems with somewhat different 
values of r. 

We averaged over an ensemble of bond configurations 
and the number of different realizations ranged from sev- 
eral hundreds to several thousands. In the simulations the 
one-cluster flip Monte Carlo algorithm was used, generat- 
ing several lO'' clusters per realization close to the critical 
point. As in earlier studies the statistical errors for each 
realization were significantly smaller than those obtained 
by averaging over the different realizations. This is the 
reason of using a relatively large number of realizations 
in the ensemble averaging. The details of the parameters 
used in the MC simulations are given in the case q = 8, 
r = 10 in Table 1. 

In the MC simulations we calculated the magnetization 
profile, defined as 



Table 1. Details of the parameters used for MC compu- 
tations (Wolff algorithm). These values are given for the 
case q = 8,r = 10. 25% of cluster flips have been discarded 
for thermalization. 



t 


Lx K 


^ of realizations 


# of cluster flips 


0.02 


640x320 


317 


10000 


0.05 


160x160 


150 


10000 




320x320 


184 


10000 




640x320 


291 


10000 


0.1 


320x320 


303 


10000 


0.15 


320x320 


303 


10000 




160x320 


84 


10000 




160x160 


56 


10000 


0.175 


160x160 


100 


10000 


0.2 


320x320 


119 


10000 


0.225 


160x160 


100 


5000 


0.25 


160x160 


187 


5000 


0.275 


160x160 


100 


5000 


0.3 


160x160 


187 


5000 




80x160 


400 


5000 


0.35 


160x160 


187 


5000 


0.4 


160x160 


187 


5000 


0.45 


160x160 


187 


5000 


0.5 


160x160 


187 


5000 



0.70 




["l(0]av = 



lX!'^'.'=l] = 



(6) 



320 



Fig. 1. Profile of the q = 6 model with free BC for two 
values of r at a distance t = 0.02 from the critical point. 
The size is 320 x 320 and average was performed over 180 
and 44 realizations of disorder for r — 8 and 10, respec- 
tively. 



where mi^k = iliSai^k-i) ~ l)/(9 ~ 1) is the local Potts 
order-parameter and the summation goes over the spins 
in the Z-th column, k = 1, 2, . . . , iiT. The brackets (. . .) 
and [. . .]av stand for thermal and ensemble averages, re- 
spectively. The absolute values are taken in order to ob- 
tain non- vanishing profiles for finite systems. The surface 
magnetization is given by [mi]av = [''7i(l)]av = \rn(V)\w 

The local magnetization, [m(Z)]av, shows a monotonic 
decrease on approach to the free surface, due to the re- 
duced coordination number close to the boundary. This is 
illustrated in Figure 1 for the random (7 = 6 model with 



dilutions r = 8 and r = 10 at the same distance t — 0.02 
from the critical temperature T^ in equation (5). The re- 
duced temperature is defined by < =| K — K^ \ jK where 
K — J /ksT. As it can be seen, the randomness tends to 
reduce order, thus the magnetization is decreasing with di- 
lution. The magnetization profile displays a plateau at the 
center of the system, the value of which defines the bulk 
magnetization, [rn^jav = ["^(i/2)]av, at the given temper- 
ature. The surface region of the profile has a characteristic 
size of ^r , which is expected to scale like to the bulk corre- 
lation length, ^ ~ t~'^ , as the critical point is approached. 
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For the random Potts model the correlation length expo- 
nent, ly, is close to 1, for all values of q [24]. 
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Fig. 2. Temperature dependence of the surface and bulk 
magnetization (g = 8, r = 10). The dashed lines are guide 
for the eyes where the extrapolated exponents of Table 3 
have been used. 



In the thermodynamic limit, L ^ cxd, as the critical 
temperature Tc in (5) is approached, the magnetization 
profile [r7T,(^)]av goes to zero as a power-law, [TO(Z)]av ~ 
tm)^ where /3(1) = /3{L) = /3i and /3(0 = /3 for ^^ < / < 
L — ^r, where f3i and /3 are the usual surface and bulk crit- 
ical exponents, respectively. The temperature dependence 
of bulk and surface magnetization is shown in Figure 2. 
To estimate the values of these critical exponents from 
simulation data one may define temperature-dependent 
effective exponents 



/3cff(0 = 



dln[m(0]a 
dlnf 



(7) 



which are approximated by using data at discrete temper- 
atures, say, t + At/2 and t — At/2. In the limit of suffi- 
ciently small At and t the effective exponents approach the 
true critical exponents, presuming that the system is large 
enough so that finite-size effects play no role (to avoid 
finite-size effects, L should be much larger than the size 
of the surface region, £^r, and the bulk correlation length. 

In the actual calculation, we approached the critical 
point by calculating /3eff(0 fo'" several temperatures, t, 
ranging from 0.05 to 0.45 with At — 0.05. As shown in Fig- 
ure 3 and 4, the effective exponents of the random q = 8 
and q = 3 Potts models approach linearly their limiting 
value. To obtain an accurate estimate for the true criti- 
cal exponents we analyzed and extrapolated the data for 
/3eff(0 using different types of correction to scaling forms. 
The most successful correction for the surface magnetiza- 
tion, written as 



TOl 



was obtained with ~ 1. 



at'^'{l + bt'^), 



(8) 
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Fig. 3. Temperature-dependent effective exponents /3cff (0 
for the surface and bulk magnetization for g = 8, r = 10. 
In the case of the bulk, the error bars are smaller than the 
symbol sizes. 



CO. 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 \- 



0.0 



° Peff(l) 



-^^~ffi$^^ 



° Peff(i/2) 



--1-. 



0.0 



0.1 



0.2 



0.3 



Fig. 4. Same as Figure 3 , for g = 3, r = 5. 



The estimates for the surface and bulk magnetization 
critical exponents are given for different dilutions in Ta- 
ble 2. The pure case values at g = 3 and 4 have also 
been computed to check the method [38]. The case g = 4 
involves known logarithmic corrections which were taken 
into account [38,39]. In the random case, as already ob- 
served in Finite Size Scaling studies by different authors 
(e.g. in Ref. [21]), due to crossover effects, the disorder am- 
plitude has a sensible influence on the exponents. At the 
optimal disorder amplitude deduced from the behaviour of 
the central charge [27] , the exponents reach their random 
fixed point values summarized in Table 3. As seen in the 
Table, both surface and bulk critical exponents depend on 
the value of g and there is a monotonic increase with in- 
creasing g. Provided that the correlation length exponent 
is close to 1 for any value of g, this observation is in accor- 
dance with previous estimates on the bulk magnetization 
scaling dimension xt = (i/v obtained in Ref. [27] at the 
random fixed point and recalled in the Table. 

At this point we are going to check the self-averaging 
properties of the local magnetization in the vicinity of the 
system critical temperature. For non-self-averaging quan- 
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Table 2. Bulk and surface exponents deduced from the 
approach to critical! ty at different disorder amphtudes. 
The variation of the exponents, outside the standard de- 
viation, is due to crossover effects. 



<? 


r 


/3i 


AP, 


P 


Ap 


3 


1 


0.541 


0.009 


0.112 


0.002 




4 


0.542 


0.010 


0.135 


0.010 




5 


0.542 


0.011 


0.1361 


0.0008 




10 


0.504 


0.020 


0.141 


0.004 


4 


1 


0.666 


0.009 


0.0831 


0.0002 




4 


0.56 


0.02 


0.1332 


0.0004 




7 


0.561 


0.022 


0.142 


0.002 




10 


0.534 


0.029 


0.146 


0.003 


6 


8 


0.581 


0.028 


0.149 


0.003 




10 


0.566 


0.018 


0.149 


0.003 


8 


10 


0.597 


0.023 


0.1513 


0.0004 



Table 3. Bulk and surface exponents deduced from the 
approach to criticality at the g-dependent optimal disor- 
der amplitude r* . The last column recalls the bulk scaling 
dimension Xb = (3 /v obtained with the same disorder am- 
plitudes in Ref. [27]. 



Pi 



P 



Xb 



3 5 0.542(10) 

4 7 0.561(22) 
6 8 0.581(28) 
8 10 0.597(23) 



0.136(1) 0.132(3) 

0.142(2) 0.139(3) 

0.149(3) 0.146(3) 

0.151(1) 0.150(3) 



tities, the reduced variance does not vanish in the ther- 
modynamic limit, indicating that fluctuations never be- 
come negligible [40] . Here we studied different moments of 
the local magnetization and determined the correspond- 
ing critical exponent, /?'•"•', defined through [m"]av ^ 

t^ [41,42,43]. For self-averaging quantities, (3^"'' is ex- 
pected to be independent of n. As seen in Table 4 the crit- 
ical exponents both for the bulk and surface magnetiza- 
tions are found to be independent of n, at least within the 
error of the numerical calculations. Thus we conclude that 
the local magnetization as the critical point is approached 
is self-averaging. This observation is in agreement with 
the expectation, that outside the critical point, where the 
system size is much larger than the correlation length, the 
central limit theorem is expected to apply, which implies 
self-averaging behaviour. At the critical point, however, 
where the above argument does not hold one may obtain 
non-self-averaging behaviour, as was observed recently by 
Olson and Young [22] for the bulk spin-spin correlation 
function. We are going to study this issue in Section 5. 



Table 4. Test of self-averaging in the off-critical behaviour 
of surface and bulk magnetization (g = 8, r = 10): The 

power law behaviours of [TO"]av and [m^Jav define only 
two exponents, for boundary and bulk behaviours, respec- 
tively. 



n 



P\ 



(n) 



P 



(n) 



0.01 0.601(21) 0.1516(4) 

1 0.597(23) 0.1513(4) 

2 0.592(21) 0.1513(5) 

3 0.587(20) 0.1513(5) 

4 0.582(21) 0.1513(5) 



4 Magnetization profile in strips at criticality 

4.1 Conformal profiles in homogeneous strips 

In a system which is geometrically constrained by the pres- 
ence of surfaces, the local order-parameter evolves from 
the surface towards the bulk behaviour and the appro- 
priate way to describe the position-dependent physical 
quantities is to use density profiles rather than bulk and 
surface observables. This is particularly important close 
to the critical point where the correlation length, which 
measures the surface region, is diverging. 

For example, in a homogeneous critical system, infinite 
in one direction, k g] — oo, -l-oo[, and confined between two 
parallel plates, which are at a large, but finite distance 
L apart, the local order parameter to(^) varies with the 
distance I £ [^,L] from one of the plates as a smooth 
function of l/L. According to the Fisher and dc Gcnnes 
scaling theory [44] : 



m(/) = L--'Fab{l/L) 



(9) 



where a and b denotes the boundary conditions at the two 
plates. In the middle of the strip, I = L/2, one recovers 
the Finite Size Scaling behaviour of the bulk magnetiza- 
tion m(L/2) ~ L"^'. In two-dimensions, conformal invari- 
ance gives further constraints on the profile: Considering 
a semi- infinite system described hy z ~ x + iy = pe , 
y > 0, with boundary conditions a and b on the positive 
and negative x axis, respectively, under the logarithmic 
transformation w(z) ~ — Inz = k + il, one obtaines the 
above strip geometry. Ordinary scaling, as in Eq. (9) then 
implies a functional form in the half-plane [45]: 



m(z) == y '"•Gab{x/p), 



(10) 



which is transformed in the strip geometry to the following 
expression 

m{w) — \w' {z)\^^''m{z) 

= [fsin(^l)]""'G„fc(cos^?/L), 

(11) 

where the scaling function Gab{cosTrl/L) depends on the 
universality class of the model and on the type of the 
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boundary conditions at the two edges of the strip. In the 
foUowing, we consider the fixed-free BC, i.e. we fix the 
spins to the state ai^k — 1 only at one boundary of the 
system, the other surface being free. As in Section 3, we 
choose periodic BC in the vertical direction, (Ti^k+i = <^i,i- 
The conformal and scaling results are strictly valid as 
ivT — > cxD, however the corrections for K ^ L are expected 
to be small. We indicate these boundary conditions by set- 
ting a = 1, & = / in equation (11). For such conformally 
invariant, non-symmetric boundary conditions, the scaling 
function has been predicted for several models [45,38]: 



Gif{cos7Tl/L) = A 



cos 



ttI 
2L 



(12) 



We mention that with the functional form in equation 
(12) one recovers the usual finite-size scaling behaviour, 
m{L) ~ L~^\ close to the free surface at criticality. The 
typical shape of the magnetization profile in a strip with 
fixed-free BC is shown in Fig. 5. 




Fig. 5. Profile with fixed-free BC in the transverse di- 
rection. In the computations, periodic BC are used in the 
longitudinal direction. 



4.2 Application to random systems 

In random systems, the conformal invariance prescriptions 
are recovered after disorder average [26,46]. The bulk crit- 
ical properties of disordered Potts models have been inves- 
tigated at criticality through conformal techniques by var- 
ious authors [23,24,25,26,27], using the longitudinal corre- 
lation function decay along periodic strips or the magneti- 
zation profiles in square-shaped systems with fixed bound- 
ary conditions along all the surfaces. Here, we consider 
the magnetization profiles of long strips with fixed-free 
boundary conditions in the transverse direction. Since a 
sufficient strip width is needed in order to apply the con- 
tinuum limit conformal results in the transverse direction, 
TM techniques are useless (the strip width is limited to 
L < 10 using the connectivity TM) and MC simulations 
are preferred (with L < 40). The parameters used in this 
work for the MC simulations are given in Table 5. We 



Table 5. Parameters for the MC simulations at criticality. 
The same parameters are used for all values of q. The num- 
ber of disorder realizations is increased at larger sizes. At 
each disorder realization, 1000 MC sweeps are discarded 
and 5000 MC sweeps are used to compute physical quan- 
tities. 



L 


K 


# of realizations 


10 


100, 200, 300, 400 and 500 


1000 


14 


100, 200, 300, 400 and 500 


1000 


18 


138, 277, 555 and 833 


1000 


24 


80 and 96 


1000 




192, 384 and 500 


4000 


30 


100 and 200 


1000 




300, 400 and 500 


4000 


40 


100 and 200 


1000 




300, 400 and 500 


4000 



mention that in spite of the large sizes used, the contin- 
uum limit is only approximately reached and perturbing 
effects will be expected close to the boundaries. 

Examples of profiles with fixed-free BC are shown in 
Figure 6 for q = 3 and 8 for different strip lengths if at a 
fixed width i = 40. The infiuence of the length K of the 
strip becomes negligible when K > 300. 



0.5 



-- ^=8 



-^■^ 



'^^l 






'<:!-^ 



0.5 1 

(/-l/2)/L 

Fig. 6. Profile with fixed-free BC for two values of q. The 
size L X K oi the strips of width i = 40 ranges between 
40 X 200 (upper profiles) and 40 x 500 (lower profiles), and 
average was performed over 1000 to 4000 realizations of 
disorder. Error bars are smaller than the symbols. 



Introducing the variable C. 



1-1/2 



in equations (11) 



and (12), one thus expects the following behaviour: 

[m(C)]av - ^(Vi)"" X [sin^]-"" x [cos7rC/2]=^i. (13) 

In order to simplify the following expressions, we introduce 
the functions /(C) = sinTr^ and g{Q = C0S7rC/2, and the 
ratio 

["^(C)]av 



i?(C,C')- 



MC% 



(14) 
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Since /(C) is symnietric with respect to the middle of the 
strip C — I' the surface dimension xi can be deduced from 
local symmetric values of the profile: 



RiCA-C) 



.9(0 



5(1-0 



[cot7rC/2]^S 



or 



Xi 



lnJ?(C,l-0 
ln[cot7rC/2] 



(15) 



(16) 



Examples of effective surface scaling dimensions, ac- 
cording to Equation (16), are shown in Figure 7 for g = 4 
with 4000 different configurations of couplings for strips 
of width L = 40 and increasing lengths from K — 100 to 
500. 
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0.5 



0.4 



0.3 



0.2 



0.1 



ii:=500- 



,t-:::s- 



- «^;"(0 



K=100 



.i=*=*^=*= 
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Fig. 7. Examples of effective surface and bulk scaling 
dimensions for q ^ A and a disorder amplitude r* = 1 . 
The size of the strip are 40 x K, with K = 100, 200, . . . , 
500. Between K = 300 and 500, the effective exponents 
remain constant up to the accuracy of the computations. 
Average was performed over 1000 to 4000 realizations of 
disorder. 



In the thermodynamic limit, these scaling dimensions 
should be unambiguously determined for any value of the 
position (^ in the transverse direction of the strip. In prac- 
tice, due to the finite size of the system and to lattice 
effects close to the boundaries ^, equation (16) defines ef- 
fective quantities xi {(, L) which do depend on the position 
C, along the strip, and also on the strip width and evolve 
towards the right limit when i — > oo. As it was already 
visible in Figure 6, between K = 300 and K — 500, the 
strip length can be considered to be long enough in or- 
der to avoid finite-size effects in the long direction. For 
these strip lengths, one can also observe a plateau in Fig- 
ure 7 where there is no significant variation of the the 



^ Close to the surfaces, (^ ^ 0, lattice effects and probable 
corrections to scaling spoil the results, while near the middle 
of the strip width, C, -^ 1/2, the precision becomes very low 
due to the proximity of the points used for the computation of 
the exponents. 



Table 6. Extrapolation in the thermodynannc limit of the 
scaling dimension of the surface magnetization, a;i(C,oo), 
measured at different values of ^. The last column presents 
our definitive determination for each value of q. 



C^o 



0.20 



0.25 



0.30 



Xl 



3 5 


0.438(1) 


0.518(1) 


0.526(2) 


4 7 


0.453(1) 


0.541(2) 


0.553(2) 


6 8 


0.478(1) 


0.567(2) 


0.577(2) 


8 10 


0.482(1) 


0.577(2) 


0.588(3) 



0.525(3) 0.523(2) 

0.552(3) 0.549(2) 

0.574(4) 0.573(3) 

0.588(5) 0.584(3) 



effective exponents which remain almost constant in the 
region ( = 0.20 — 0.35, and whose extrapolated values 
should be consistent. The values computed in the plateau 
region are thus studied as L increases and extrapolation 
is made towards the thermodynamic limit L — > oo. The 
dependence of the effective exponents on the strip width 
is shown in Figure 8 for C = 0.25. A linear extrapolation 
leads to the estimations of a;i(C,oo), collected in Table 6. 
The results exhibit a good stability relative to the position 
C and allow a final determination of the scaling dimension 
Xl, given in the last column of the table. On the other 
hand, we observe important corrections to scaling in the 
boundary behaviour, since the surface exponent falls down 
rapidly as C ^ 0. This effect is the possible origin of the 
discrepancy with the boundary scaling dimension found 
in Ref. [19] for q = 8: xi ~ 0.47, which is recovered here 
in the limit (^ —f 0. 
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Fig. 8. Effective surface exponent at C = 0.25 plotted 
against L~^. The intercept corresponds to the extrapo- 
lated value in the thermodynamic limit: a;i(0.25, oo). 



For the bulk exponent, using the quantity 

R{Cl/2)^[f{C)]-^n^9iOr\ (17) 
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we can form a combination where xi cancels, leading to: 



Xb 



{lni?(C,l/2)xln[V2g(l-C)] 
-lni?(l-C,l/2)xln[V2.9(C)]} 



{ln[/(C)] 



X In 



9(C) 



9(1-0 



(18) 



This expression involves three points of the profile and is 
thus subject to stronger numerical fluctuations than the 
surface scaling dimension, as shown in Figure 7. Unfortu- 
nately, compared to previous precise determinations, no 
accurate extrapolation can be made here. One can how- 
ever check the conformal expression (13). For that pur- 
pose, it is necessary to extrapolate the rescaled profiles 
["T'l(C)] X -^^S obtained at finite sizes, towards the ther- 
modynamic limit. Considering, as we did before, that the 
longer strips are large enough to be unaffected by finite- 
size effects in the long direction, extrapolation to infinite 
width L only will be done. For each position ( in the 
transverse direction, we plot the corresponding local mag- 
netization as computed for different strip widths and then 
extrapolate in the limit L ^ oo. Examples of linear and 
quadratic least square fits are shown in Figure 9. The lat- 
ter one has been preferred. 
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0.02 



0.04 0.06 

1/L 
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Fig. 9. Extrapolation to infinite width of the profiles. 
The solid lines corresponds to q = 8 and the dashed lines 
to q = 3. Two fitting curves are presented (for ( — 0.25, 
g = 3), a least square linear fit (dotted line) and quadratic 
(solid line). The quadratic fit is more accurate. 



It leads to an extrapolated profile, which is compared 
to the conformal expression in Figure 10 for q — 8. In 
the conformal formula, the scaling dimensions of Tables 3 
and 6 are entered and the amplitude is the only free pa- 
rameter. In spite of the small strip widths considered, the 
extrapolation procedure introduced in Fig. 9 is very accu- 
rate, since the agreement between extrapolated data (full 
circles) and Eq. 13 (solid line) is fairly satisfactory. This 
is a strong evidence which supports the validity of the 
conformal expression for the profile. 



g 0.8 






I — I 0.4 




Fig. 10. Rescaled profiles at different strip widths for 
q = 8 (open symbols) and extrapolation to L — ^ oo (full 
circles). The solid line is the conformal expression where 
the amplitude is the only free parameter. The agreement 
between the extrapolated data and the conformal profile 
is very good. 

5 Multifractal behaviour at criticality 

In this section, we study the possible multifractal behaviour 
at the critical point and consider different moments for the 
magnetization profile, [m"(/)]av , and that of the corre- 
lation function [G'^(r)]av • The characteristic exponents, 

x]^"' and xi'' , in the bulk and at the surface, respectively, 
are expected to vary with n for multifractal behaviour. 
Indeed, it was first Ludwig [8] who predicted multifractal- 
ity in the bulk correlation function of the random bond 
Potts model by conformal perturbative methods (see also 
the results by Lewis in Refs. [47,48]), which was confirmed 
recently by Olson and Young [22] by MC simulations in 
the square geometry. Here we rather work in the strip ge- 
ometry and calculate both the bulk and surface scaling 
exponents. 

The calculations about the moments of the critical pro- 
file are parallel with that in the previous section and the 
scaling dimensions are extracted from the expected func- 
tional form: 



["^"(C)] 



1/n 
av 



(") 



A{TT/Lfb X [sinTrC]"'^'' x [cos7rC/2]"=i 



(19) 

which is analogous to (13) for the average behaviour, i.e. 
for n — I. We note that the typical behaviour corresponds 
to n = 0,i.e. exp[ln?n(C)]av. 

In the actual calculation we have considered the q = 
8 model on strip-like samples with K = 500 and L = 
10, 20, 30 and 40 and the average is performed over 5000 
realizations. Using the method of the previous section first 
we deduced from equation (19) the effective, size and po- 



sition dependent surface scaling dimensions. 



(") 



(C,i), 



which are then extrapolated to the thermodynamic limit, 
L —f oo. The extrapolation procedure is demonstrated in 
Figure 1 1 , whereas the extrapolated data are presented in 
Table 7. As one can see in this Table the critical point 
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Table 7. Extrapolation in the thermodynamic hmit of the Furstenberg method [50] : 
scahng dimension associated to the moments of the surface 

magnetization and measured at different values of ( (q — ^ 

8). The last row presents our definitive determination. ^o(i) = lim — In 



surface magnetization of the random bond Potts model 
shows multifractal behaviour: The scaling dimensions of 

(n) 

the different moments of the surface magnetization, x^ , 
are monotonously decreasing with n. Note that for n = 2 
the surface scaling dimension is very close to the pure (and 
random) Ising value of xi — 1/2. 
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Fig. 11. Effective surface dimensions of the n -order mo- 
ments of the magnetization profile, evaluated at C = 0.25 
for q = 8 and n = 0, 1, 2, 3. Dashed lines denote results 
of quadratic extrapolation. 



m^oo rn 



Rt. 



\vo} 



(20) 









C 






n 


0.20 


0.25 


0.30 


0.35 


4"' 



1 
2 
3 


0.601(7) 
0.585(5) 
0.496(3) 
0.387(11) 


0.600(7) 
0.582(6) 
0.496(4) 
0.384(14) 


0.602(9) 
0.582(8) 
0.500(6) 
0.384(18) 


0.599(11) 
0.579(11) 
0.498(8) 
0.381(25) 


0.600(9) 

0.582(8) 
0.498(5) 
0.384(17) 



where T, is the transfer matrix and 



is a unit initial 



vector. The free energy density is thus given by [/o(i)]av = 
— L~^ Ao{L) . For a specific disorder realization, the spin- 
spin correlation function along the strip 



G,(fc) 



q{S, 



0'lA<^l.k + l. 



-1 



q-l 



(21) 



follows from the application of products of transfer matri- 
ces on the ground state eigenvector associated to Aq (for 
details see, e.g. Ref. [27]). 

We will now assume that conformal covariancc can be 
applied to the order parameter correlation function and 
its moments. In the infinite complex plane z = x + iy 
the correlation function exhibits the usual algebraic decay 



nl/n _ 



^^" where 



at the critical point [G^(p)]a4" = const x p 

p = I zi — Z2 I . Under the logarithmic transformation w = 

2^ In z — k + U , one gets the usual exponential decay 



a;(")fc 



along the strip [G"(fc)]av" — const x exp ( — ^ 

Fig. 12 for an illustration). The scaling dimension a;^ 
thus be deduced from an exponential fit. 



(n) 



see 



can 




Next we turn to study the multifractal behaviour of 
the bulk magnetization. As mentioned in the previous sec- 
tion, from the magnetization profiles in finite strips one 
cannot extract precise estimates for the scaling dimen- 
sion xj^^' . Therefore we used a different technique, based 
on the Blote and Nightingale connectivity transfer ma- 
trix [49]. Since transfer operators in the time direction 
do not commute in disordered systems, the free energy 
density is defined by the leading Lyapunov exponent. For 
an infinitely long strip of width L with periodic boundary 
conditions, the leading Lyapunov exponent is given by the 



Fig. 12. Spin-spin correlation function on a long strip 
with periodic BC in the transverse direction (the upper 
left corner also corresponds to the upper right corner). 
This Figure is a 3D plot where the value of the correlation 
function between the origin (upper left corner) and any 
lattice site {k, I) is plotted on a vertical axis at position 
{k,l). 



This method was used in Ref. [27] for the average cor- 
relation function, i.e. for n = 1. In this work, we calcu- 
late the higher moments, as well as the typical behaviour, 
which corresponds to n = 0. 
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For each strip size {L = 2 — 8), we considered sys- 
tems of length ~ 10® and an average is performed over 
80 X 10'^ disorder configurations. For a given strip the ef- 
fective, size-dependent exponents follow from a linear fit 
in a semi-log plot, as exemplified in Figure 13. It is clearly 
seen that the scaling dimensions are different for the dif- 
ferent moments of the spin-spin correlation function. The 
effective exponents are then extrapolated as L — > oo and 
the estimated values are presented in Table 8 for different 
values oi q (at r*{q)). 



10" R, 
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40 

2%k/L 



60 



80 



Fig. 13. Moments of the spin-spin correlation function 
(semi-log scale) ioi q = 8, L ~ 7. 



Table 8. Scaling dimensions of the moments of the spin- 
spin correlation function computed at r = r*. 







4 


n) 




n 


q = 3 


q = A 


q = Q 


g = 8 



1 
2 
3 

4 
5 


0.154(1) 
0.132(1) 
0.116(1) 
0.104(1) 
0.095(2) 
0.088(2) 


0.177(1) 
0.138(1) 
0.114(1) 
0.097(1) 
0.087(2) 
0.079(2) 


0.207(1) 
0.146(1) 
0.114(1) 
0.094(2) 
0.081(2) 
0.072(2) 


0.234(1) 
0.150(1) 
0.112(1) 
0.091(2) 
0.077(2) 
0.068(2) 



Again the critical point bulk magnetization shows mul- 
tiscaling behaviour, for all q > 3 the scaling dimensions, 
Xjj" , are monotonously decreasing with n. One can make 
a comparison with pcrturbative results: 



and A ^ 11/12 - 41n2 and B = ^(33 - 29V37r/3) are 
numerical factors. The first-order term in Eq. (22) is due 
to Ludwig [8] and the second-order term to Lewis [47]. 
This result is known to be valid close to q = 2 and for the 
lower moments. The comparison with our data is made in 
Table 9. One can observe a good agreement with the re- 
sult of Lewis for q = 3, especially for the lower moments. 
For q — A, the perturbation expansion probably breaks 
down already at n = 3, since xj^ is found to be larger 

than Xjj with the formula of Lewis. It is remarkable that 
the second-order term in the perturbation expansion van- 
ishes at n = 2, which could explain the stability of the 
exponent associated to the second moment, with respect 
to variations of q, has already been noticed by Olson and 
Young [22]. 



Table 9. Scaling dimensions of the moments of the spin- 
spin correlation function computed at r = r*. The nu- 
merical estimates on the first column are compared to the 
results of the 1st and 2nd order perturbation theory in 
Eq. (22). 





x^" for the 3-state Potts model 


n 


num. 


1st order 


2nd order 



2 
3 


0.154(1) 
0.116(1) 
0.104(1) 


0.158 
0.108 
0.083 


0.157 
0.118 
0.110 




x]^ for the 4-state Potts model 


n 


num. 


1st order 


2nd order 



2 
3 


0.177(1) 
0.114(1) 
0.097(1) 


0.188 
0.063 
0.000 


0.181 
0.120 
0.167 



We close this section by presenting different moments 
of the magnetization profile, [rn^ {l)\JJ^ , as extrapolated 
towards L ^> oo. In Figure 14 two moments (n = and 
n = 2) of the scaled profiles are plotted for the g = 8 
model, where the open symbols represents the finite-width 
data. The extrapolated profile, denoted by full circles is 
in perfect agreement with the conformal results in equa- 
tion (19), where the scaling dimensions x^ and x]^ are 
taken from Tables 8 and 10, respectively, and the only 
fitting parameter is the amplitude in equation (19). We 
consider this agreement as a strong evidence in favour of 
the validity of the conformal expression for the averaged 
moments of the order parameter profile. 



^(") _ ^p 



9 



=xr^("-i) 



+ (A + B(n-2))e2 + 0(e3) 



(22) 

where a;^ is the bulk exponent in the pure model, e is a 
small expansion parameter related to the deviation of the 
pure model's central charge to that of the Ising model 
(e = 2/15 and 1/3 for q — "i and 4, respectively [48]), 



6 Discussion 

In this paper the critical behaviour of the two-dimensional 
random bond Potts model is studied by MC simulations, 
transfer matrix techniques and conformal methods. New 
features of our present work are the following, i) For the 
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Fig. 14. Rescaled n'^'-order profiles at different strip 
widths for q = 8 (open symbols) and extrapolation to 
L ^ oo (full dots). The solid line is the conformal expres- 
sion with only the amplitude as a free parameter. 



first time we have investigated the surface critical be- 
haviour of the model and determined the surface magne- 
tization critical exponent, /3i, from approach to criticality 
In addition we got estimates for the corresponding bulk 
exponent, /3. ii) We have studied the critical point magne- 
tization profiles in strip-like geometries and deduced the 
scaling dimension xi = Pi/v from the predictions of con- 
formal invariance. in) We have presented numerical ev- 
idence for the multifractal behaviour of the surface and 
bulk magnetizations at the critical point and the scal- 



ing dimensions of the averaged moments. 



J") 



and Xi 



(n) 



are calculated, iv) Finally, we have demonstrated that the 
different moments of the critical magnetization profiles, as 
well as the correlation function obey conformal invariance. 
The critical exponents and the scaling dimensions of 
the average quantities are continuously varying functions 
of q, however an exponent relation 4/3 = j3i is approxi- 
mately satisfied, as can be observed in Table 3. The anoma- 
lous dimensions of the relevant scaling fields in equations 
(1) and (2) can be estimated using the scaling relations: 
Vt = 1/j^ = Xb/13 = xi/fii, y/i = 2 - Xb, and yh^ ^l-xi. 
Their values are presented in Table 10. While yt remains 
close to 1 (but V satisfies the limit v > 2/d [51]) for all 
values of q, the anomalous dimensions related to the mag- 
netic field vary with q. 



Table 10. Anomalous dimensions of the relevant bulk and 
surface scaling fields. 



Xb/P 



Xl/Pi 



Xb 



■ X\ 



3 5 


0.971(29) 


0.965(6) 


1.868(3) 


0.477(2) 


4 7 


0.979(35) 


0.979(13) 


1.862(3) 


0.451(2) 


6 8 


0.980(40) 


0.986(18) 


1.854(3) 


0.427(3) 


8 10 


0.993(26) 


0.978(18) 


1.849(3) 


0.416(3) 
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